Inelastic scattering effects and the Hall resistance in a 4-probe ring 
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Phase randomizing processes in mesoscopic systems can be described in a phenomenological way 
within the Landauer-Biittiker formalism by attaching extra voltage probes to the sample. In this 
paper, it is shown that a perturbation treatment of this idea allows for the incorporation of such 
effects without the need of giving up the efficiency of recursive techniques commonly used for calcu- 
lating the transmission coefficients. The technique is applied to a 4-probe ring, where a Hall effect 
can be observed that originates from quantum interference rather than a Lorentz force acting on the 
electrons. The influence of inelastic scattering on both the Hall resistance and the Aharonov-Bohm 
oscillations in the longitudinal resistance are examined. 



I. INTRODUCTION 



II. MODELING INELASTIC EFFECTS 
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Quantum coherence plays a prominent role in the 
transport properties of mesoscopic systems; interference 
between different electron trajectories can lead to inter- 
esting effects like weak localization 1 and Aharonov-Bohm 
(AB) oscillations. Inelastic scattering events will destroy 
phase coherence and as such interference effects will be 
smeared out. Modelling inelastic scattering proves to be 
difficult because it is in general a many-body problem. 
However, a proposal has been put forward some years ago 
for incorporating these effects phenomenologically within 
the Landauer-Biittiker (LB) formalism (which is a single 
particle theory), by attaching extra voltage probes to the 
sample-. Nevertheless, only few papera&i have made use 
of it because the widely used standard recursive Green's 
function method is unable to calculate the transmission 
coefficients between the extra voltage probes. As such, 
one should in principle invert the complete Hamiltonian 
to solve for all transmission coefficients which is of course 
computationally very inefficient for large systems. In this 
paper it will be shown how to treat the regime of weak 
inelastic scattering with a perturbation approach to the 
original voltage probe model, still keeping the computa- 
tional effort the same as needed for the standard recursive 
Green's function method. 



Our method is applied to a ring setup with 4 probes 
and a flux piercing through like depicted in Fig. ^ I 11 
such a setup a Hall voltage can be observed which is 
completely due to interference of electrons travelling in 
opposite directions along the ring, and which in princi- 
ple does not rely on a Lorentz force acting on the elec- 
trons 5 ' 6 . In the current paper, numerical results for a 
realistic two-dimensional (2D) ring will be shown, taking 
into account inelastic scattering. The observed Hall effect 
disappears when decreasing the phase coherence length, 
showing again that the effect is solely due to quantum 
interference. 



For calculating resistances in our system, we will use 
the LB formalism which relates resistances to transmis- 
sion probabilities between the leads. In the formalism, 
the leads are thought to be connected to large reservoirs 
with a well-defined chemical potential and temperature. 
The currents through the leads and the voltages on the 
reservoirs are related through (at temperature T = 0): 



(i) 



where p, q label the leads/reservoirs and T pq is the trans- 
mission probability from lead q to lead p. Although this 
formalism is only valid when phase coherence is present in 
the sample and the leads, phase breaking events must be 
taking place in the reservoirs in order for them to have a 
well-defined equilibrium distribution with a certain chem- 
ical potential. It is this insight that was used by Biittiker 
to arrive at the idea that an extra lead connected to a 
reservoir can introduce a phase breaking events. Indeed, 
if the current through such a lead is fixed to be zero (this 
is what we call a voltage probe), then for every electron 
that enters the lead and is absorbed by the reservoir, an- 
other one has to come out. But since equilibration is 
taking place in the reservoir, the electron coming out is 
not coherent with the one going in. In this way, it is 
possible to model inelastic effects in a phenomenological 
way. 

In our calculations, we will attach a one-dimensional 
(ID) lead on every site of our tight binding model in order 
to have a homogeneous distribution of inelastic scatter- 
ing centers. These leads can be thought to extend in a 
direction perpendicular to the 2D sample (see Fig.0. An 
adequate choice of the potential in the ID leads makes 
it possible that the influence of these leads can be de- 
scribed by just adding a complex onsite energy — ir) to 
every site in the lattice^. The parameter r\ is controlled 
by the hopping parameter between the ID leads and the 
sample, and is related to the inelastic scattering time as 
■q = H/2t v . However, the approach goes further than 
just adding an imaginary potential; in order to conserve 
the total current in the system, one has to solve Eq. (JTJ 
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so that the current through the extra ID leads is zero. 
Before giving expressions for the currents through the 
leads, we will introduce some notation: the ID voltage 
probes attached to the system to simulate phase ran- 
domizing processes are labelled by a Greek index, while 
the real physical leads connected to the sample (referred 
to as contacts from now on) will be labelled by letters 
m, n, . . .. Now, by putting I a = for all a, and solv- 
ing for the voltages V a on the voltage probes in terms 
of the voltage differences V n — V m between the contacts, 
one can obtain an expression for the currents /„ through 
the contacts similar to Eq. QJ, but now with effective 
transmission probabilities incorporating the effect of the 
voltage probes: 



In = 



2e 2 
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where the effective transmission probabilities can be writ- 
ten as follows: 
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with S-y = J2i T-yi+J^g^ T-yS- This expression has a clear 
physical interpretation: the first term describes transmis- 
sion from contact m to contact n without any inelastic 
process, the next term incorporates a single scattering 
event, then double scattering and so on. 

A standard method commonly employed (because of 
its efficiency) for calculating the transmission coefficients 
is the recursive Green's function technique (reviewed in 
Ref. 0). However, this technique can give only access to 
the transmission coefficients between the contacts T nm . 
For calculating the effective transmission coefficients in 
Eq. ©, one would also need the transmission probabil- 
ities between the contacts and the voltage probes (T na 
and T am ), and between the voltage probes themselves 
{Tap)- These are unavailable with the standard tech- 
nique, because one would need Green's functions con- 
necting points at the edges of the sample (where the 
contacts are attached) with points inside the sample, 
whereas the technique only gives us Green's functions 
between the left and right edge. 

However, one could neglect all terms in Eq. J2J involv- 
ing two and more scattering events (by putting T a p = 0) , 
and keep only the direct transmission term together with 
the single scattering term. This approximation is valid 
when the phase coherence length of the sample is larger 
than the system size. Although the standard technique 
still cannot give access to T na needed for the second order 
term, a recursive technique has been developed recently 
that allows to calculate these with the same numerical 
effort needed for the standard recursive Green's function 
method. Details of this technique can be found in Ref. 8. 
In that way we are able to treat the effective transmission 
coefficients to second order in the phase randomization 
processes. 




FIG. 1: Schematical view of the 4-probe ring setup. Inelastic 
scatterers will be modelled by attaching a ID lead at every 
scattering center. In real calculations, every site of the tight 
binding model is attached to a voltage probe in order to have 
homogeneous scattering. 



It should be noted that by neglecting higher order 
terms in Eq. (j3J) , one will break current conservation rules 
whenever a magnetic flux is present in the system; but in 
the weak scattering regime, this error will be small. In 
all calculations we present, the sum of currents flowing 
through leads 1 to 4 divided by the sum of their absolute 
values was smaller than 10~ 4 . 



III. RESULTS FOR A 4-PROBE RING 

We will consider a 4-probe ring in a two-dimensional 
electron gas, like shown schematically on Fig.^ All leads 
are arranged to the left and right of the ring for compu- 
tational convenience. The Hamiltonian of the system can 
be written as: 



H = — (p - eA) 2 + Mmp + Vo, 
Zm* 



(4) 



where m* and e are respectively the effective mass and 
charge of the electron, A is the vector potential created 
by the flux through the ring, V c is the confinement po- 
tential defining the ring, and Vi mp is the elastic impurity 
potential. By discretizing this Hamiltonian on a square 
lattice with lattice parameter a, one obtains a tight bind- 
ing model: 

H = y^y^e mra |w,n)(TO, n\ + 

n tci 

(d,nl m i 71 + 1)( TO ; n \ + *m.rJ m + 1' n ) { m i n \ + ^' C -) 

with (m, n) labelling the sites on the lattice, and e mn the 
on-site energies. The hopping parameters are given by: 



_ te -ie/h [ A.d\ 



(5) 



with the integral evaluated along the hopping path and 
t = h 2 /2m*a 2 . 

The ring parameters were chosen in accordance to typ- 
ical rings fabricated experimentally in a 2DEG at the 
interface of an GaAs-AlGaAs heterostructure. The den- 
sity of the 2DEG was chosen to be n s — 4 x 10 11 cm~ 2 , 
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FIG. 2: The Hall resistance R H = 1/2 (#12,34 - #34,12) for 
different values of the phase coherence length. 
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FIG. 3: Longitudinal resistance of a 4-probe ring for different 
values of the phase coherence length. 



corresponding to a Fermi wavelength of 40 nm. The 
ring has a mean radius of 0.6 /itm, and the width of 
the arms is 200 nm, so that 10 channels are available 
for conduction. The mobility of the electron gas is 
fx = 5 x 10 5 cm 2 V -1 s _1 , giving a mean free path of 
5.2 fim. 

This translates into the following tight binding param- 
eters. The lattice parameter was chosen to be a = 6.7 nm, 
so that Xp — 6 a and the Fermi energy Ej = 1.1 1. 
The mean radius of the ring is 89 a, and the width 29 a. 
Elastic scattering was modelled with on-site energies 
distributed randomly in the interval [— 0.127£, 0.127i], 
which gives a mean free path of l m — 780 a (estimated in 
Born approximation) . 

For this ring, we will calculate both the longitudinal re- 
sistance i?i2,i2 and the transverse resistances i?i2,34 and 
#34,12 ( see Fig-DJ, and the influence of inelastic processes 
on the results will be made clear. We use the common no- 
tation Rki,mn = {V m — V n ) / Ik for a measurement where 
current is supplied through contacts fc and I, and the volt- 
age difference V m — V n is measured, fixing I m = I n = 0. 
In terms of transmission coefficients, these resistances are 
given bjai: 

nkl.mn = Y) ' W 

where D is a 3 x 3 subdeterminant of the matrix A relat- 
ing the currents through the four contacts to their volt- 
ages [I m = A mn V n , c.f. Eq. D is independent of the 
indices klmn. When inelastic processes are included, the 
transmission coefficients in these expressions have to be 
replaced with the effective probabilities given in Eq. 

Let's look at the resistance #12,34- Suppose an elec- 
tron enters the ring through lead 1. It can reach lead 3 
by different paths: there is a direct path between lead 
1 and lead 3, but there is also a trajectory going like 
1 — > 4 — > 2 — > 3. These two trajectories will interfere 
with each other, and fix the voltage on lead 3 (neglecting 
paths circling the ring more than once) . The same can be 



applied to lead 4: again there is a direct and indirect path 
interfering and setting the voltage on lead 4. It is clear 
that because of symmetry, the voltage on lead 3 and lead 
4 will be the same when no flux is present through the 
ring. However, when a flux pierces through the ring, time 
reversal symmetry is broken and the phase difference be- 
tween both paths going from 1 to 3 is different from that 
for the paths going from 1 to 4, so that a voltage differ- 
ence V3 — V4 will develop. For a one-dimensional ring, this 
effect is described in some detail in Ref. 0, and numeri- 
cal calculations for a simple model are shown in Ref. 6. 
The effect is purely based on interference; a Hall voltage 
will also be present when the magnetic field is limited to 
the inside of the ring, and no flux is going through the 
arms of the ring. Like in the original Aharonov-Bohm 
proposal, the electron can feel the vector potential, and 
is not subjected to a real Lorentz force. 

The Hall resistivity is the anti-symmetric part of the 
resistivity tensor: pn = l/2(p xy — p yx ). In our case 
this comes down to calculating the Hall resistance Rh = 
1/2 (#12,34 — #34,12), which is the transverse resistance 
part that is anti-symmetric with respect to the magnetic 
nu>si. Results are shown in Fig. [21 The Hall resistance 
varies periodically with the magnetic flux through the 
ring, the period being the fundamental flux quantum "Jo- 
in the figures, only one fundamental period is shown. It 
should be mentioned again that there is no flux through 
the ring arms, so the flux is fully contained within the 
ring. The interesting effect is that for nonzero flux, a 
Hall resistance can be measured, which is not due to a 
Lorentz force, but which originates from interference be- 
tween different trajectories along the ring like explained 
previously. 

Phase randomizing processes are included in the way 
explained in the previous section. The phase coherence 
length of the sample is given by L v = yj Dt v , where 
D = l/2vFl m is the diffusion constant, and t v = h/2rj is 
the phase relaxation time. The parameter r\ corresponds 
to the hopping parameter between the attached voltage 
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FIG. 4: Hall resistance Rh of the 4-probe ring when a mag- 
netic field is applied across the whole sample. $ is defined as 
the flux through the mean radius of the ring. The effect of 
inelastic scattering is shown in the inset, for the first 2 oscilla- 
tions in the resistance: L v = oo (solid line) and L v = 7.8/im 
(dashed). 



can be observed experimentally quite easily nowadays, 
this means that it is also feasible to measure the Hall 
effect in the 4-probe ring. 

However, in an experimental setup it is difficult to con- 
fine the magnetic flux to a region inside the ring, and 
therefore the magnetic field is applied across the whole 
sample. In Fig. 21 we show calculation results for this 
case. The Hall resistance Rh is not anymore strictly 
periodic with respect to the flux. Nevertheless quasi- 
periodic oscillations are visible resulting from quantum 
interference, whose amplitude will decrease when intro- 
ducing inelastic scattering (see the inset of Fig. [3J. A 
beating pattern can be observed in the oscillations, which 
is a result of having several channels open for conduction 
in the ring arms; trajectories for different channels sur- 
round slightly different areas and thus fluxes. Compared 
to Fig. |3 the Hall resistance also gets a non-zero offset: 
this contribution is a result of the Lorentz force acting 
on the electrons when there is a flux present in the ring 
arms. 



probes and the sample. For our system, we varied the 
phase coherence length between 6/im and infinity, which 
is clearly in the range of validity of our approximation 
(system size > L v ). When the phase coherence length is 
reduced (dashed and dotted lines in the Fig. [5J , one can 
see that the amplitude of the oscillations decreases. This 
was to be expected because the Hall effect we observe is 
completely due to quantum interference. 

Turning to Fig. |3J the normal AB oscillations are ob- 
served in the longitudinal resistance. Furthermore, one 
sees that the influence of inelastic scattering on the lon- 
gitudinal resistance is twofold. Firstly, the mean value of 
the resistance will increase by decreasing the phase coher- 
ence length because we are introducing extra scatterers. 
Second, the amplitude of the AB oscillations decreases 
with decreasing coherence length, because these oscilla- 
tions are the result of interference effects. 

When comparing the magnitude of the Hall resistance 
in Fig. [21 with the longitudinal resistance shown in Fig. [31 
we see that the oscillations have the same order of mag- 
nitude. Since normal Aharonov-Bohm (AB) oscillations 



IV. CONCLUSION 

Incorporating the effect of phase randomizing pro- 
cesses in a phenomenological way can be done with the 
attachment of extra voltage probes to sample, an idea 
originally proposed by Biittiker 2 . In this paper, we have 
described a method for treating inelastic effects based on 
this idea, but such that a recursive technique for calculat- 
ing the transmission coefficients can still be used. Doing 
so, one is able to treat the problem in a numerically very 
efficient way, which was not possible within the original 
proposal. The approach however consists of neglecting 
multiple inelastic scattering and is therefore only valid in 
a regime where the phase coherence length is larger than 
the system size. 

The method was applied to an experimentally realiz- 
able ring with four attached contacts, and a Hall effect 
was observed which is due to quantum interference rather 
than an implicit Lorentz force acting on the electrons. 
This Hall effect disappears with decreasing the phase co- 
herence length. 
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